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Abstract 



We present an approximate NNLO evaluation of the QCD form factor resumming large logarithmic pertur- 
bative contributions in semi-inclusive heavy flavour decays. 



1 Introduction 



In this note we present an approximate next-to-next-to-leading order (NNLO) evaluation of the QCD form 
factor resumming threshold logarithmic contributions in semi-inclusive heavy flavour decays. 

In semi-inclusive processes, final giuon radiation is strongly inhibited in the phase space regions where the 
observed final state obtains its maximum energy, therefore opening the way to soft and coUinear singularities. 
The perturbative calculation of the differential cross section is plagued, in that limit, by large logarithms. In 
order to improve the reliability of the perturbative calculation, these large logarithms need to be resummed. 

Let us consider the rate of semi-inclusive heavy flavour decays. The Mellin iV-moments of the rate contain 
double logarithmic contributions and have an expansion of the form: 
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where 



L = In TV. 



(3) 



Rn (as) is a remainder function, which does not contain large logarithms and has a perturbative expansion of 
the form: 



Rn {as) = a^s ^ 



•AT • 
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The running coupling is evaluated at a general renormalization scale jj, ^ Q, where Q is the hard scale: 

as = as (/i^) . (5) 
The logarithmic terms have an exponential structure, so one can write |^: 

1 



-Tn (as) = G (as) /n {as) 
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where the form factor (a^) reads: 
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exp [ci2Q!5 L"^ + ciias L + C23a| + C22a| + C2 lO?! L + cs^ag ~\ ] 
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Note that the exponent contains only the first term as of the double-logarithmic series (a^L^)". The 
advantage of the exponentiation is therefore that we can predict In /at reliably for asL <C 1, that is, for a larger 
region than asL"^ <C 1, where the perturbative expansion of /at holds. In fact, the other terms of (asL^)" come 
purely from the expansion of the exponential function, as can be seen in formula (0). 



The prefactor in (^ is the coefficient function, having an expansion in powers of as'- 

00 
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The double sum in the exponent is usually organized as a series of functions, which resum "strips" in the (n, k) 
plane: 
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exp [L gi iPoasL) + g2 (PoasL) + as 53 iPoasL) + ag .94 (PoasL) 
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The functions gi (A) have a power-series expansion: 

00 

ff.:(A) = 5]5»A". 
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where A = fi^asL. They are all homogeneous functions: gi (0) = 0. This property insures the normalization 
of the form factor: fN=i — 1- The resumming at LO, NLO and NNLO is referred, respectively, to series in 
L(asL)", (asi)" and asiasL)^. 

We have computed the function (73(A), which is necessary for the NNLO resumming. Not all the quantities 
determining (/3(A) are known exactly. A3 is only known numerically from a fit of the three-loop Altarelli-Parisi 
(AP) splitting function to the known moments. B2 is unknown and we approximate it with the infrared-regular 
part of the two-loop AP splitting function. The coefficient functions and the remainder functions for radiative 
and semileptonic B decays have been computed to NLO in ref. A complete NNLO computation of the 
distributions requires also the knowledge of the two-loop coefficient function, which at present is unknown for 
any distribution. After the resummation of the threshold logarithms in the TV-space, we have returned to the 
form factor in the a;-space, by inverting fN{as) with analytic and numerical procedures. 



2 QCD form factor at NNLO 



Let us briefly describe the derivation of the NNLO form factor in the iV-space. By NNLO accuracy, we mean 
the resummation of all the infrared logarithms up to and including 
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The general expression of the form factor is: 



log/jv (as) = / dz 
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Let us remind one difference between annihilation processes, as Drell-Yan, and other processes, like the 
present one or DIS. In the annihilation processes, the initial partons reduce their momenta by irradiation, 
before actually annihilating; when t = Q^/s — > 1, only the emission of soft partons is allowed. Likewise, in 
the other processes, as f.i. DIS, and in the limit x = Q^/2Q ■ p ^ 1, only soft emission is allowed, before the 
scattering; however, after the scattering, the parton fragments with the only kinematical constraint of having a 
low virtuality and coUinear emission is no longer forbidden. 
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The functions A (ag) , D {as) and B (as) have a perturbative expansion in powers of as' 



ri=l 
oo 

D{as) = ^Dna'^ = Dias + D2al + --- ■ 

n=l 
oo 

B {as) = ^" = Bias + ^2^1 + • ■ • . (13) 



Let us observe, that in general the transverse momentum rule is a guess, as it has not been proven to such 
an accuracy. The universality of B2, in general, is a debated problem. If we neglect the variation of the coupling 
with the scale (frozen coupling) , we find logarithmic terms of the form: 

Alasi^ A2alL^, Aaa^L^,--- D^asL, Dzali, • • • B^asL, B2a|L, • • • . (14) 

Then, to NNLO accuracy, one needs the first three terms in the expansion of A {as) and the first two terms of 
the functions D {as) and B {as) ■ 



The three-loop coupling, according to the definition given in the PDG ji), reads: 



I /3i loglogA^VA^ /3i Mog'log>VA'-loglogA^VA'-l P2 1 

/?olog/.VA2 (31 logVVA^ 131 logWA2 logVVA^ ' 



The asymptotic expansion of the coupling is basically an expansion in inverse powers of log/x^/A^. The first 
three coefficients of the /3-function, defined as 



da„ 



d log ^2 



-/3oa| - Pial - P2a% , (16) 



read: 



f3o = ^- - = — — — ^ = 0.87535 - 0.05305 ni.. 

Pi ^ - ^^^"f - ^^--^ ^ - ^ 0.64592 - 0.08021 n., 

^ 2857 - 5033/9 n. + 325/27n^ ^ ^^^^^^^ _ ^^^^^^^^^^ ^ ^ ^^^^^^ ^2 . (,7) 
1287r'' 

Let us note that /3o and /?i are renormalization-scheme independent, while P2 is not and we have given its value 
m the MS scheme [||. 

Integrating the /3-function equation (|l6|) on both sides, one obtains: 



as (fc') = as (Q") - /3o4 (Q') log 7^ - /3i4 (Q^) log — - /?2a| (Q^) bg + (iterations). (18) 



P 1,2 1,2 

Substituting the above expression into (|l2|), one sees that a /3o-insertion corresponds to the additional factor 
as L, that of (3i to a| L and that of P2 to ag L : 

Po ■■ as L, Pi : a| L, P2 : a| L. (19) 
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Therefore, in the terms containing NNLO coefficients, the coupling can be replaced with the one-loop one and 
in the NLO terms the couphng can be replaced with the two-loop one, so that one has: 



A{as) = Aias,3L + A2al2L + ^30^s.il'^ 
D{as) = Dias,2L + D2alj^^; 
B{as) = BiasaL + B2al,y^. 



(20) 



Furthermore, the j3\ term in A2a| 2L '^^^ tie neglected, as it is a N'^LO contribution. After replacing eqs. ( [20| ) into 
eq. (12), one performs a straightforward integration over fc^, the gluon transverse momentum. The integration 
over 0, the longitudinal gluon momentum, is easily done using the approximation Q]: 



l-z 



(21) 



This approximation misses the term proportional to >liC2, where C,2 = ""^/S, which can be obtained in the 
following way. Using the large- iV approximation derived in one obtains, in the n variable: 
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Therefore: 



[ 1 1 ^2 52 

[nnlo] = — 



2 9 (log 



[lo] 



(23) 



We have used the variable n ~ N/Nq instead of N, where A^o = e^^'^ — 0.561459 .... Within this alternative 
representation, the terms proportional to and to jj^ disappear. This scheme is probably more accurate 
as Feynman diagram computation directly in A^-space brings factors containing F (1 — e) with _D = 4 — 2e the 
space-time dimension. 

The advantage of the variable A^ is that the total rate is directly reproduced by setting A^ = 1 , while in the 
variable n it is given by /„=i/jVo- These two variables differ by terms of higher order in je- 



The lowest-order term, computed within the approximation (^), reads: 

A\ ' s s s s 

P°] = ~T7:r log ^ log log ^ + log s log log s - 2 log - log log - 
2po L n n 



(24) 
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One then obtains: 



log /at (as ) = - — - [log s/n^ log log s/n'^ + log s log log s - 2 log s/n log log s/ n] 
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The next step is to write the above result as a function of the three-loop coupling. The inverse of ([l5| ) reads: 



logs 



^,^L_,A,„„,w-(f-|)o., + o(4), 



(26) 



where s is the square of the hard scale in unit of the QCD scale, 

After the replacement (p6|), terms of higher order with respect to NNLO a re g enerated, which must be discarded. 
One finally makes the expansion in up to second order to obtain eq. ([40|). 



(27) 



The quantities A\ and Ai are known analytically Q 

Cf 



Ai = — = 0.424413, 

TT 



= 0.42095- 0.03753 nF, 



(28) 



where Ca = = 3, Tr = 1/2 and ni? = 3 is the number of active quark flavours in h decay. The value given 
for Ai in the MS scheme for the coupling constant. At present, only a numerical estimate of A3 is available 



A3 = 0.59413 - 0.09272 np - 0.00040 n\. 
The soft quantities D\ and D2 |0] are known analitically: 



(29) 



TT 
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' ^2 
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TT Tr n_F 



-0.59826 - 0.07157nF. 
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Numerically, C (3) = 1.20206. The constant Di is renormalization-scheme independent, while D2 is not and we 
have given its value in the MS scheme. The latter quantity turns out to be the most important one to determine 
the NNLO effects. 



The coUinear quantity Bi is known analitically. 



^ --^ = -0.31831. 

4 TT 



(31) 



The two-loop quantity B2 is unknown and we approximate it with the infrared-regular part in the two-loop 
Altarelli-Parisi splitting function Pqq (z). In general, the latter is naturally decomposed in the soft limit as: 



P,, (^) 



1 - z 



- K (z; as) + K (as) S {I ~ z) . 
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The functions K (z; as) and K (as) have an expansion in powers of as: 



(32) 



K (as) — Kias + K2as + ■ ■ ■ , 
K{z;as) = Ki (z) as + K2 (z) a^ + 



(33) 



Since the splitting function has vanishing first moment: 



it holds that 



K (as) = / K{z;as)dz, 
Jo 



Ki (z) dz and K2 



In leading order, it holds Bi — Ki. Our approximation then is: 



K2 (z) dz. 
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where, in the MS scheme |10 
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2.1 Results 

The functions gi and 172 have the following expressions |l^, |l^ : 
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-f - [(1 - 2A) log (1 - 2A) - 2 (1 - A) log (1 - A)] 
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+ ^ [log(l - 2A) - 21og(l - A)] + ^ [log(l - 2A) - log(l - A)] + 
^Po PO 

-T^ [log'(l - 2A) - 21og2(l - A) + 21og(l - 2A) - 41og(l - A)] + 

+ ^ log(l - 2A) + |l log(l - A) + A [log (1 - 2A) - 2 log (1 - A)] log ^. (39) 
2po Po 2po 



6 



Our result for the NNLO function reads: 
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(40) 



Arbitrary constants have been added to the function 93 in order to make it homogenous. The quantity 7^ = 
0.577216 ... is the Euler constant and ( (n) is the Ricman zcta function, 



00 -. 



k=l 



(41) 



C (2) — TT^/G — 1.64493. The functions 172 and gs depend on the renormalization scale /U, while gi does not. 
The expansion up to order reads: 

log/jv = —AiasL'^ - DiasL- BiasL + 

-Ia,(3o a| - (]-A2 + PoDi + IpoBi) a| - ( D2 + B2 + Ia,C2Po] a^L + 



--Ai/3o'a|L* - (^A2/?o + ^A^Pi + ^Dif}'^ + r^B^Paj Oi%L' + 
- + \a^C2PI + 2i?2/3o + B2/30 + + \biP^ a%L^ + O (4) • 

Let us note that /?2 appears only at order O (a^) . 



(42) 
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The functions gi become singular when 



Since 

A = Poas (a*^) 
this means that a singularity in iV-space occurs when: 

N ■ 



exp 



A' 



(43) 



(44) 



(45) 



where A is the QCD scale. Let us observe that, in general, this singularity signals non-perturbative effects but 
its precise position is completely unphysical, as we can move it with a change of renormalization scale. 



2.2 Renormalization-scale dependence 

In this section we consider renormalization-scale dependence. In principle, such scale fj, should not appear in 
the cross sections, as it does not correspond to any fundamental constant or kinematical scale in the problem. 
The completely resummed perturbative expansion of an observable^ is indeed formally independent on /i. In 
practise, truncated perturbative expansions exhibit a residual scale dependence, because of neglected higher 
orders. 

We start with the form factor as a function of 0^5 (Q^) and we derive its expression as a function of as (/i^) 
and fJ-^/Q^. 



Since 



a (Q) = a {fi) +ca^ (^) + c' (fi) + O (a^) , 
where a = Poas and (from eq. (p^)) 



c-log^ 



log + — 



one has, where now A — A(/i): 



(46) 



(47) 



Lgi[\ + acX + a^c'A] = L gi [A] + cX^g'i [A] 
(?2 [a + acA + a^c'A] — g2 [A] + acXg'2 [A] 



+ ac'X'g[ [X]+a-c' A Vi' [A] 



The additional terms in the functions gi, to (partially) compensate for the scale change 



(48) 

/i^ therefore read: 
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1. d 
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2 M2 
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2 M 



Pi 



log- 



A^ (A^ g[ [A]) log^ ^ + (Ia^s; [A] + Xg', [A]) lo^ 



A .9^ [A]log- 



(49) 



^We do not mean here the resummation of towers of logarithmic contributions, but the resummation of the whole series in ag. 
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Therefore, the leading function gi does not depend expHcitely on the renormahzation scale, analogously to the 
coefhcient of the leading term in the expansion of inclusive observables (such as the ag/Tr term in the total e~^e~ 
hadronic cross-section). The higher order functions gi^i instead explicitly depend on fi. We want to show the 
effect on log[7V] of the variations in eq. (^). Therefore, in fig. Qand in fig. ||, we show the dependence of log[iV] 
on /i at NLO and at NNLO, keeping as fixed at 0.21; we notice a NNLO improvement at high values of N. At 
lower values of N (N < 10), there is an appreciable improvement only by going at much higher energies (f.i. 
as = 0.1). In general, by varying also as with the scale /i around rrih, there is not a significant improvement at 
NNLO. 



3 Analytic inverse Mellin transform 



The iV-moments are physical quantities, but in practice a measure of the moments for large N is difficult. It is 
therefore convenient to perform the inverse transform back to momentum space. 

Let us derive an analytical expression of the inverse Mellin transform at NNLO. We start by introducing 
the inverse Mellin transform of fN{as)/N: 



1 /-O+Joo ^j^ 

G{as;x) = — —x~^fN{as,L) 



(50) 



The inverse Mellin transform of /at (as) is the logarithmic derivative of G{as, L; x) 

I{x)^-x^G{x). 



(51) 



After the substitutions x — e " and u = Ns, we have: 



G{as;x) 







du 


27ri 


J C — ioo 


U 




/•C'+ioo 


du 


2TTi 


'c'~ioo 


u 



fN{as,L) 

OO 

e- exp L gi (PoasL) + ^ ag~^ gn {fioasL) 



n=2 



(52) 



We Taylor expand the exponent with respect to L around i = Z = In 1/s |1J 
ln(— Inx). Note that I — ln(l ~ x) when x ^ 1. 



We have, at NNLO: 



G{as;x) = — 
Ztt I 



C +200 



^ gU+-Fo(0+-Fi(01n"+^^'2(0 ln= u+. 



C'-ioo 



where, at scale and at the same NNLO: 

Fa{l) = 'gi(/3oasO + 52(/3oa50 + "s53(/3oas0> 
F2{1) - 2l3oasg[{,l3oasl)+l3lallg'^{l3oasl). 



In the x- variable, I 



(53) 



By keeping in the exponent terms up to NLO, while expanding the NNLO terms, we obtain: 



G{as\x) 



1 



1 + F/" ^(/) \nu+-F2{l) In u + 



(54) 



9 



where: 



9\ iPoOisl) +Po(Xsl g'l {Poasl) 
Ff^il) = l3oasg'2{l3oasl). 



(55) 



By using the result: 



27riJc ^ rfi^iVLfc r(l - i^f ^ 



)' 



(56) 



where F is the Euler Gamma function, we obtain, after the integration: 



G{as;x) = 



1 + Ff^ V(l - i^i^"^) + ^ F,il) (V^(l - Fl"^) - ^'{1 - ^)) 



r(i - Ff^) 

Qlgi(l3oasl)+g2(0aasl)+as gaWoasO 

T{l-gi{f3oasl)-poaslg[{f3oasl)) 



[l + f3oas g2{Poasl) 1^(1 - giiPoasl) - Poas I g'liPoasl)) 



+- F2{1) - giiPoasl) - fioaslg'iiPoasl)) - - gi{Poasl) - Po^s I g[{Poasl))) 



(57) 



with il>{x) = dlogr{x)/dx, the digamma ftmction. 

At the end, we find the expHcit analytic formula for the inverse Mellin transform at NNLO: 



fix) 



-x^G{as;x) 



= —X- 



d C qI gi(0oasl)+g'2{fiaasl)+as gzifioasl) 

{ r(l - giiPoasl) - Poas I g'liPoasl)) 



[1 + Poas g'2{Poasl) - giiPoasl) - Poas I g[{Poasl)) 



+- F2{1) (V'(l - giiPoasl) - Poaslg'iiPoasl)) - ^'(1 - giiPoasl) - Poas I g[{Poasl))) 



(58) 



3.1 Discussion and conclusions 



There are two physical effects involved in the inverse Mellin transform: 



1. exact longitudinal momentum conservation; the terms which enforce longitudinal momentum conservation 
are formally subleading in the infrared logarithm counting; 

2. infrared pole in the coupling; the inverse transform involves an integration over all moments and arbitrarily 

large values of |A^| enter. This implies that a prescription for the infrared pole in the coupling has to be 
given. In general, the transformation mixes all the momentum scales in the problem. 



One can study one problem at a time; for example, one can study the problem 1. only, by considering the 
frozen coupling approximation for the distributions in A''-space. In one-loop, this is equivalent to consider the 
QED case. 
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The frozen coupling approximation means neglecting the variation of 0:5 with the scale. Let us start from 
formula (|l|) at NNLO: 



rQ^{l-z) ^^2 

log fN (as) - / dz— / — [Aias + A2al + Asal + ...] 




Bias + B2al + ... + Diag + -D2a| + ...+} 

/•I z^-i -If 1 ^ 

J dz—^—— I {Axas + Aaal + Ay.al) In ^— + {Bi + Di) as + {B2 + D2) a|| 



(59) 



After integration in z, at the lowest order in A = /3o'^2 In A^, we have: 



„ = I _|i _ 52 _ i|I£'| X (62) 



/3o Po Pa 
B2 D2 A2J1 
/3o /3o /5o 

We have two ways, analytical and numerical, to compute the inverse Mellin transform of /at (as): 



1 



2TTi 



f{as;x)^—l dNx-"" fnias^L). (63) 



C—ioo 



In the case of frozen coupling, the gi are linear in A and therefore the numerical path does not include the 
Landau pole; the numerical integration is therefore exact. We have compared the inverse Mellin transform 
calculated numerically, directly from the definition (|63|), with the analytical result (|58|). We find a very good 
agreement, up to energies around the charm scale, provided that we go to NNLO. The exact longitudinal 
momentum conservation does not seem to spoil the reliability of the perturbative series. 

In the real world, we cannot use linearized Qi, but we have to use the gi computed in section |2| therefore, 
problems due to infrared poles come into play. 

Let us first make an important observation: the degree of singularity of the functions gi for A 1/2, and 
therefore also of the form factor, increases with the order of the function, i.e. with i [Q. 

At LO, gi has at most a logarithmic singularity of the form 

<?i: (l-2A)log(l-2A); (64) 

at NLO, (72 has at most a singularity of the form 

<?2 : W (1 - 2A) . (65) 

The prefactor 1 — 2A in front of the logarithm, vanishing for A = 1/2, is now absent, but the singularity is still 
logarithmic. At NNLO, 53 has at most a singularity of the form 

log' (1 - 2A) 
1-2A ■ 

The latter is basically a pole singularity, no more a logarithmic singularity anymore, i.e. a much stronger 
singularity. Note that the above term is proportional to Ai0l and so it is scheme- independent. 
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Another observation concerns the size of the two-loop soft term. For three active flavours, the size of the 
two-loop correction with respect to the one-loop one is rather large, 



^as « 2as « 40%, (67) 

since as {ttib) — 0.21. The inclusion of D2 is therefore important. The size of the two-loop term compared to 
the one-loop term is expected on general ground to be of order: 

^ « ^ = 2.25, (68) 

since in two-loop order gluons start to be radiated by gluons istead of quarks, the former having a larger colour 
charge Ca instead of Cp- 



Let us now consider the inverse Mellin transform (|63|), releasing the frozen coupling approximation; in 
formula ( [l^ ) the coupling runs over the whole integration range. 

The form factor in A^-space is computed in such a way that infrared-pole effects appear in a sharp way for 
N > Nc^ where it acquires an unphysical (imaginary) part. In other terms, the numerical distribution is not 
real for any value of N because of the integration over the Landau pole. An exact numerical evaluation of the 
inverse transform then requires a prescription for the pole. An alternative strategy is to give a prescription for 
the infrared pole directly in A-space, in such a way that the form factors are well-defined for any A^. In general, 
this results into a softening of the form factor for large A^. It is then not necessary to give a prescription for the 
pole in the inverse transform. 

We have used the minimal prescription (mp) |l5[| , over two different paths (with the same results); the first 
path was made by two straight lines parallel to the negative real axes, closed by a half-circle centered arond 
the origin and crossing the positive axes between the origin and the first Landau pole; the second path was 
composed by two lines almost vertical, meeting on the positive real axes between the origin and the first Landau 
pole. The precise crossing point is irrelevant, as far as it is before the Landau pole. 

The inverse Mellin transform can be also derived analytically, as seen in section |[ We want to compare the 
analytic and numerical results, as we did before in the frozen coupling approximation. In this case, however, 
the Landau pole contribution plays an important role and we need to reach regions of higher energy {as — 0.1), 
where the first Landau pole (see eq. (^3|)) occurs at larger values of A^ (TV = e^/'^>^"°"' ~ 3000 - 4000). In these 
regions, the perturbative resumming keeps under control the infrared divergences and cancels the oscillatory 
behaviour. We compare the analytic (NLO and NNLO) and numerical (NNLO) plots in fig. || (at 0:5 — 0.1). 
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Figure 2: We show the dependence of log [AT] on the scale fi at NNLO, with ~ ms (solid line), fi ~ 2m b 
(dashed line) and /i ~ 1713/2 (dot-dashed line), at fixed asiQ"^) = 0.21 and 5 < N < 25. 
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Figure 3: mp numerical distribution (solid line), NLO (dashed line) and NNLO (dot-dashed line) analj^ic 
distributions of /(a;), at = 0.1. 
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